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The rotation of the polarization vector for light propagating perpendicular to an external constant 
external magnetic field B, is calculated in quantum vacuum, where it leads to different photon 
eigenmodes of the magnetized photon self-energy tensor for polarizations along and orthogonal to 
B (Cotton-Mouton effect in QED vacuum). Its analogies and differences with Faraday effect are 
discussed and both phenomena are calculated for a relativistic electron gas at low densities, by 
starting from the low energy limit of the photon self-energy eigenvalues in presence of B. In the 
QO . Cotton-Mouton case the polarization vector describes an ellipse whose axes vary periodically from 

zero to a maximum value. By assuming an effective electron density of order 10 3 cm -3 the quantum 
relativistic eigenvalues lead to a rotation of the polarization plane compatible with some of the limit 
fSJ i values reported by PVLAS experiments. Other consequences, which are interesting for astrophysics, 

are also discussed. 
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I. INTRODUCTION 



For light propagating in magnetized quantum vacuum, in a direction perpendicular to the magnetic field B, and 
polarized along a plane forming a nonzero angle (3 with B, the polarization vector rotates describing a curve which 
\ is a sort of ellipse whose axes oscillate in time, making the polarization plane to rotate. This effect is due to the 
birefringence of magnetized quantum vacuum (for a medium it is known as Cotton-Mouton effect), which exhibits 
different dispersion laws for light polarized along and perpendicular to B. In the present letter we obtain the frequency 
of the ellipse axes oscillations by starting from the difference of the photon self-energy tensor eigenvalues, which we 
' calculate by starting from the low frequency limit of the expressions obtained by Shabad [l[ for the photon self-energy. 
t-H , We find a similar phenomenon for a photon propagating in a relativistic electron gas, and from the low density limit of 
the corresponding eigenvalues [|J we obtain the oscillation frequency of the ellipse axes, leading to the Cotton-Mouton 
effect. By a Lorentz boost parallel to B it results that the rotation occurs for any photon non-parallel propagation 
00 ' direction, whose polarization is also non-parallel to B. 

We will discuss in the present paper both cases of propagation in quantum vacuum and in a dilute electron gas, 
which may arise from the ionization of neutral atoms or molecules. But, based on general arguments, even in the 
case that such background charge is negligible, a dilute neutral gas (for instance, that remaining in the ultra-high 
vacuum achieved in the laboratories, or in astrophysics, the components of some gaseous nebulae) would exhibit the 
birefringence effect, if atoms and molecules are under the action of an axially symmetric external field B. Their 
interactions with photons would differ for polarizations parallel and perpendicular to B. As a result, it is expected 
■ to have different frequencies for these two directions, leading to a Cotton-Mouton rotation. This effect bears some 
analogy to the Faraday rotation, due also to birefringence, which occurs in a magnetized medium non-invariant under 
charge conjugation for photon propagation parallel to B. Thus, in QED vacuum Faraday effect does not exist, but 
it occurs if a charge background is present. For a charged relativistic electron -positron gas its dispersion modes and 
polarizations were found in |2|, and some additional features are discussed below in Section III. However, the Cotton- 
Mouton rotation occurs in vacuum as well as in any other media with no restriction about the charge conjugation 
property. Our results are interesting both in connection to PVLAS experiments and in astroparticle physics. 

In magnetized QED vacuum the spatial symmetry is explicitly broken by the field B. Electrons and positrons 
(observable and virtual) move in bound states characterized by discrete Landau quantum numbers 0,1,2.. on the 
plane orthogonal to B (the quantum version of the classical circular motion), and move freely along B. By using the 
appropriate vectors characterizing the reduced symmetry properties and from gauge invariance,(see Q) one can 
write the general tensor structure of the photon self-energy or polarization operator , which is different from that 
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of isotropic vacuum (B = 0). The polarization operator can be written as a linear combination of three basic tensors, 
which are even functions of the electromagnetic field tensor F^ v . 



II. LOW ENERGY PHOTON EIGENMODES AND POLARIZATION ROTATION IN VACUUM 

To understand what follows, we recall briefly some properties of the photon eigenmodes in a magnetic field B 
The diagonalization of the photon self-energy tensor leads to the equations Ii^ v Cv^ — K^C/p having three non 
vanishing eigenvalues and three eigenvectors for i = 1, 2, 3 (One additional eigenvector is the photon four momentum 
vector k v whose eigenvalue is — 0). The first three eigenvectors satisfy the four dimensional transversality 
condition = 0), and are C* = k 2 F 2 x k x - k^(kF 2 k), C% = F* x k x , C 3 = F^ x k x (Here F^ = d^A u - d v A„ 

is the electromagnetic field tensor and F* u — ^t^vpnF 9 *" its dual.). By considering a^\x) as the electromagnetic four 

vector describing these eigenmodes, its electric and magnetic fields are = —g^Q^ — ^ a o \ — V x aM . The 
polarization vectors are given in detail in @ , [l[ . 

The eigenvalues and eigenvectors characterize the eigenmodes of wave propagation in the magnetized vacuum. In 
what follows we consider sometimes the eigenvectors normalized to unity (which is valid whenever the magnetic field 

B is not taken as zero), and denote them by a£ . We must recall here that for propagation perpendicular to the field, 
the second mode has polarizations proportional to = b\\(k 2 — Lu 2 ),h^ = —{b± x by]. We use the set of unit 
vectors bi ; 2 , (or in general bj_) orthogonal to B, and b3 = b|| along it, whereas the third mode has polarizations 

e^ 3 ) = — [bj_ x byjwjhj 3 = — bxfcii. The first mode is purely longitudinal e^ 1 ' = — hj_uj and non-physical in quantum 
vacuum. In a charged medium it is responsible of longitudinal waves. For propagation along B, the modes 1,3 are 
transverse with polarizations proportional to = —hj_ui and e^ 3 ) = — [bj_ x bu]u>, respectively, whereas the second 

is pure longitudinal = b||(fc| — lu 2 ) and consequence it has meaning only in a charged medium. The vectors kj^ 
and k|| are the components of k across and along B. The previous formulae refer to the reference frame which is at 
rest or moving parallel to B. 

These results indicate that the dispersion equations and polarizations of light in a magnetic field are different along 
different directions in space. Let us consider a plane polarized wave E = fe luJt propagating perpendicular to B and 
whose polarization vector f is in a direction forming a nonzero angle [3 with B. It can be decomposed in two waves 
having polarizations f 2 = tcosfl, f 3 — fsin{3 along and perpendicular to B.Thus, after entering in the magnetized 
quantum vacuum, each one of these components obey the dispersion law dictated by the corresponding eigenmode, and 
propagate with different polarization and frequency; birefringence is produced. The frequencies are LU2 = lu + Aa>2 and 
co>3 — lu + A<x>3, where Aw2,3 are the contributions from the corresponding eigenvalues of the polarization operator 
«2,3 respectively, obtained after solving the dispersion equations k 2 — «W — 0. At a point of space we can write the 
projections of the electromagnetic wave E as Ei = f2Cos(uj + Aa^i and £3 = fzcos(iu + Ao^i. The polarization 
vector starts to rotate, since we have here the well known problem of the resultant curve from orthogonal oscillations 
whose frequencies differ in a small fraction of the main frequency. 

We will give below a quantitative estimate of the rotation frequency. In the low frequency limit, under the additional 
condition k\ « 2eBhc, the scalars and have been approximated (see Appendix) by expanding the general 
expression of the photon self-energy tensor obtained in [l[ and taking the first order in B 2 . The expansion is made 

k 2 

perpendicular to B we have z\ = —lu 2 .) 



also in terms of the relativistic invariant variables z\ = k 2 — lu 2 , kj_, respectively as (notice that for propagation 



« (2) = ^r(- 2 + k)> (i) 
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These quantities contain the contribution of all Landau levels of virtual electron-positron pairs. The frequency e 
varies in the present approximation as B 2 . On the light cone k\ + z\ = 0, (p}, © agree with earlier values reported 
by Dittrich et al [T3] in calculating the vacuum birefringence in a strong magnetic field. 

The estimated ratio of frequencies e/to = (k,^ — K^ 3 ')/k±cu for lu ~ 10 15 rad/s, is c/lu ~ 10~ 21 . 

To find the curve described by the polarization vector we shall write lu' = lu + Alu2 for simplicity, and by eliminating 
lu' from the expressions for E2 , E3 one gets 



f 2 Ej - 2f2f 3 E2E 3 coset + f 2 E 2 = flflsm 2 et 



(3) 



FIG. 1: The ellipse described by E2, E3 is drawn here for three values of /3, along the time axis. We have exaggerated the ratio 
e/u to ~ 10~ 2 , to make the rotation of the polarization plane more evident. 

The curve representing © is an ellipse with oscillating axes, expressed in terms of E2, E3, and /2, fz- For et = 0, 7r, 
it degenerates in two straight lines E 3 — ±tan/3£?2- For e = 7r /2, 3tt/2 we have respectively two ellipses. The ellipse 
starts with its minor axis being zero and its major axis along the initial plane of polarization. The major axis oscillates 
and decreases, while the minor axis oscillates increasing, up to a configuration in which the major axis is zero and 
the minor axis is maximum and passes along the other straight line. Starting from this configuration, the situation is 
reversed. (See Fig. [T]). This effect makes the rotation difficult to be observed in some intervals of time, when one of 
the semi-axes is very close to zero. We will see below that if some adequate density of electrons exists in the medium, 
the rotation of the polarization plane is produced at a frequency e much larger than in the quantum vacuum, allowing 
it to be observable at magnetic fields and frequencies of the same order as those of PVLAS experiments. 

III. THE CASE OF A MEDIUM 

The case of photons propagating in a relativistic electron-positron (e^) medium at temperature /3 _1 7^ and non- 
zero particle density (chemical potential /i 7^ 0) is of especial interest |2[. The one-loop diagram describing the process 
accounting for the photon self-energy interaction contains, in addition to the virtual creation and annihilation of the 
pair, the process of absorption and subsequent emission of one photon by the electron and/or positron. Concerning 
the average particle densities, for /i ^ but = 0, the average positron density is zero. This is usual in most 
laboratory conditions. 

Also, even in the case of high vacuum conditions of the experiment 0], we expect that a density of order 10 8 
molecules remain. This is far from being pure quantum vacuum, but a medium which is not invariant under charge 
conjugation. We assume the existence of a low density electron gas may due to ionization of a fraction of these 
molecules which might be produced both by the action of the laser beam and the magnetic field intensity. In such 
case, their equilibrium with ions might be described by some chemical potential /1 7^ . The ions also contribute 
to the self-energy tensor, but due to their large mass, their contribution is much smaller than that of the electron 
gas. In such a medium, an additional longitudinal electromagnetic wave component is also possible, saturating the 
three spatial degrees of freedom. Two other tensors, odd in the chemical potential fj, and anti-symmetric in the 
tensor indices p, v contribute to n p „ through new basic scalars named V, R. The resulting eigenmodes are, in general, 
polarized elliptically, and in some cases are combinations of transverse and longitudinal waves • 

For propagation along the external field B we name K2 — S m and = T m ± R, where S m ,T mi R are scalars 
depending on temperature and density as well as on and B. The second mode is a pure longitudinal wave, whose 
electric polarization vector e^ 2 ) = (uj 2 — k 2 )h^ 7^ since kj_ — but lu 2 — k 2 does not necessarily vanish in the medium. 
The transverse modes are 

C'^ = R(Cl± l Cl) (4) 
and describe circularly polarized waves in the plane orthogonal to B having different eigenvalues. Up to a normalizing 
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factor, their electric polarization vectors can be written respectively as e' 1 ' 3 ) = (b^ ± ib% ), corresponding to the 
eigenvalues T m ± R, leading to the well-known phenomenon of Faraday effect. Let us consider two electromagnetic 
waves whose polarization vectors are proportional respectively to e^ 1,3 ', (h*- 1,3 -* =kx e' 1,3 )). 

If we also assume that T mi R oj, one can write approximately uj^i = \/k 2 — T m =p e', where e' = i?/2\/k 2 — T m . 
Thus, a superposition of both modes having equal amplitudes leads to the following wave 

E = ARe[e 1 e i ( k -*- ul - t '> +e 3 e l(k - x -^] (5) 
= A[e 1 e~ u ' t +e 3 e"' t ]i?ee j(k - x -" t) 
= A(b x cos e't + b 2 sin e'^Ree 1 ^-^] 

which shows that the polarization of the propagating wave rotates counterclockwise with frequency e', describing 
a circumference. We recall that as the system has rotational symmetry around B, the direction of the orthogonal 
eigenvectors b 1 ' 3 may be taken arbitrarily. 

Concerning the propagation perpendicular to the field, we have the eigenvalues «2 = S m , Kx,3 = [{P — T m ) ± 
v/(P -T m ) 2 -4R 2 ]/2 which describe two waves whose polarization vector rotates in the plane orthogonal to B. 
However, for very low charge densities, the quantities P and R are negligibly small as compared with T m for perpen- 
dicular propagation. This means that the eigenvalues Ki 3 can be approximated as k 3 , respectively, the 3-rd mode 
being plane polarized orthogonal to B. Thus, a wave entering in the magnetized medium propagating orthogonal to 
B, with polarization vector forming some angle ir/2 > (3 > with B, has components along the polarizations of modes 
2, 3. This means that the polarization vector would describe an ellipse with oscillating axes, like the curve mentioned 
previously in the vacuum case. Below it is shown that its frequency of rotation is approximately a linear function of 
the electron density, N e . 

For comparison with PVLAS results, we take frequencies 10 ~ 1, 8 x 10 15 rad/s, and magnetic field B ~ 5.5 x 10 4 G, 
and we may simplify the expressions given in the Appendix for T m and F. Let us call the density of particles by 
N e = ^2 n (eB /h 2 c 2 ) J dp^n e , where n e is the electron Fermi distribution function. We must stress here that n e stands 
for the low temperature limit of the sum of electrons plus positrons densities n e + n p , (n p is negligible small at low 
temperatures), since the Cotton- Mouton rotation exists in any case, even in a neutral electron-positron gas, where 
n e = n p as different from the Faraday rotation, which depends on the net charge n e — n p . We may approximate (|10[) 
and (T2"j) given in the Appendix as 

T m = -^A 3 m 2 C 4 V e (6) 
S m = --A 3 m 2 c 4 V e , (7) 

since (fTUl) and (TT2"|) under the conditions of propagation perpendicular to B, k\ <C m 2 c 4 and near the light cone 
(k± ~ u>) are approximately independent of the photon momentum and energy. Here A is the electron Compton 
wavelength. The dependence on B is contained in the term N e (see below) . We assume that to have charge neutrality 
one must add to ([5J [7]) a similar term in which the electron density is replaced by the ion density and the electron 
mass by the ion mass. Such term, however, is al least 10 -3 times smaller than expressions (J5]),([T]), and may be ignored 
in the calculations. 

For the previous values of B and lj, we have finally Alu — (S m — T m )jk^ ~ ^A 3 m 2 c 4 iV e /fc^ ~ 5, 6 x 10 _34 iV e . For 
N e ~ 10 3 , e ~ 10~ 31 er<7 = 10 -4 rad/s, leading to a period T ~ 10 4 s. The angle rotated in a length of one meter is 
of order 10 -13 rad. 

Eqs. ([6|),([7]) are valid for any value of the electron density and energy. Their dependence on the magnetic field 
is contained in the density of particles term: for constant chemical potential and zero temperature, the dependence 
on B comes from the degeneracy factor eB and from the energy eigenvalues E n , especially for nonzero temperature. 
However, another dependence might come from the chemical potential fj,, which may be -B-dependent. 

Thus, for having a Cotton-Mouton rotation comparable with the limits imposed by the last results of PVLAS 
experiments [l3| , if on the average, and in the process of the experiment, the system behaves as if having a number 
of free electrons of order 0.001 percent of the remnant molecules. This is enough to give a figure 10 3 times the pure 
quantum vacuum effect. 

The Fermi-Dirac degenerate distribution might not be suitable to describe n e in some cases, for instance, at nonzero 
temperature and very low densities. Thus, although the experimental PVLAS conditions are made at temperatures 
near zero, the observed phenomenon occurs far from equilibrium. For a fixed magnetic field we may assume that on 
the average it behaves equivalently to an equilibrium system having a constant chemical potential and an effective 
temperature of several Kelvin degrees, and use the Boltzmann distribution n e = ^ a„ J dp^,e~^ E ~^^ («„ = 2 — <5q„), 
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as an approximation. For E we take the non-relativistic limit by expanding it as E = mc 2 + p 2 /2m + eBhn/mc. We 
call A — 13 , where fj! = mc 2 — p, and gets easily 

„ eB , eBHB 

N e = A- — — coth (8) 

Xonc mc 



where \u = hj yj2itm/ (3 is the De Broglie thermal wavelength. For for // constant and X < 1, by taking cothX ~ 
(l/X + X/3), we conclude that N e , and in consequence e grows with B as a + bB 2 , where a = m/X]jh 2 ,b = eh(3/mc are 
temperature dependent parameters. For T of order of few Kelvin degrees, the dominant term is a which is independent 
of B. But the chemical potential might be dependent on B, which is more realistic for the ionization process. Thus, 
for the dilute gas under the influences of both the laser and magnetic field it is expected for T m , S m a dependence on 
B more complex than in the quantum vacuum case (i.e. not proportional to B 2 ) 

We must observe that formulae (|10ll2p in the Appendix may be applied even in the case of the so-called hot vacuum, 
a system of electron-positron pairs in equilibrium with photons at high very high temperature (/3 _1 > mc 2 ~ 10 9 °K), 
Such a system have zero net charge, and is invariant under charge conjugation. Such system (which might exist in 
neutron stars magnetospheres) does not exhibit the Faraday rotation but it shows the Cotton-Mouton rotation of the 
polarization plane. Even more, as hot vacuum implies fi = 0, and the Boltzmann distribution is a good approximation 
for it [l2| one can use the expressions below (|6l7p and ((HJ) by taking A — 2e~ l3mc . It is seen from ([8]) that by increasing 
temperature, the density of particles, and in consequence, the frequency induced by the magnetic field, decreases. 

In concluding, we want to emphasize that the Cotton-Mouton rotation of the polarization plane for photons prop- 
agating orthogonal to B, although bears some analogy to Faraday effect, differs from it in two main facts: Faraday 
rotation is circular (which is related to the axial symmetry determined by the preferred direction of magnetic field 
B, whereas the Cotton Mouton rotation is a more complex curve; Faraday rotation occurs only if the system is not 
invariant under charge conjugation, whereas Cotton-Mouton rotation does not require non-invariance under charge 
conjugation. Thus, it occurs in QED vacuum, in a medium containing an electron-ion system, as well as in a neutral 
molecular or electron-positron gas. Even if the remnant ionized molecules density would be negligible small, the 
medium would show birefringence properties, with the consequent rotation of the polarization plane. The Cotton- 
Mouton rotation is expected to act on the polarization plane of radiation propagating across intergalactic media. 

The authors thank the Abdus Salam ICTP support under OEA-ICTP Net-35 Grant. They thank M. Chaichian, 
A. Tureanu and A. Zepeda for several comments and discussions, to A. Gonzalez, B.Rodriguez and S. Villalba for 
valuable suggestions, and especially to A.E. Shabad for several important remarks. 



IV. APPENDIX: 



a) We write the explicit expressions for S and T in quantum field theory in which it is taken h = c = 1, and in 
statistics [2]. In the first case, by calling z 2 = k\, the scalars n 2 ^ = S,T are given by 



where 



= — / dr d V [p 2 , 3 C]e M -£e-^ 
n Jo Ji 



eB 1-7] 2 1,1 

C = —r—B,N = —J-, M = --z 2 M' - - Zl N 
smhei? 4 2 2 



1 1 ± n . sinh cBrry sinh eBrri 
rt = -,M' = - 



eB sinh eBr 



p 2 = cosheBrN - ^z 2 Q 



P3 = -\z 2 CM' - ^ Zl Q 



_ sinh eBrrr cosh eBrri + 

V '■ 

sinh eBr 



(9) 
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One can expand ([9]) in powers of z%, Z2- To first order one gets 

K 2,3(^l, z 2) — z l\—^ ) Zl =z 2 =0 + Z 2(—K )z 1= z 2 =0, 

since K2, 3(0,0) = 0. By taking the expansion in powers of B 2 also to first order, one gets eqs. ([1]), ((2|). 

b) We shall now write from [2|,|9| the explicit expressions for the scalars R, T m and S m for cu 2 h 2 <C 2eBhc, in the 
case of a charged medium, i.e. , in relativistic quantum statistics. We have 

Tm = -^f £ F Sl f f[l - G]K + n p - 1) (10) 



n.ir 



J 'l 



where G = [(ff + iV^/IQI 2 , where = zi + 2eBhc{n + n')), K = z x + 2eBhc(n' - n), \Q\ 2 = K 2 - 4uj 2 E 2 , and 
the electron and positron Fermi-Dirac distribution functions. The quantity —1 subtracted to them stands 
for the quantum field limit /i — > 0, /3~ 4 = 0, whose contribution to order B 2 was discussed previously. The quantity 
E q — m 2 ^ + c 2 p\ + 2eBhcn is the electron-positron energy eigenvalue. The quantities F^,, N = 4jV 1 * 1 ^, are 

defined in as well as F^^) and G^ n ,, used below, in terms of the Laguerre functions of the variable k 2 _/2eB. 
They have the property F ( ^) (0) = S n , n >-! ±<5„_i,„, F^,(0) = +£„,„/, G^ n ,(0) = 2(<W +S 0n ),N^ n , (0) = 

0- 

Now we give the corresponding expression for R. We have R = <J^F, where 



eBhc 
' 4tt 2 



(11) 



The term i* 1 arises due to the electron-positron charge asymmetry. It vanishes in the limit fi — > 0, but is nonzero in 
absence of positrons, whenever the electron density n e ^ 0. 

The expression for S m , which in the present case contributes to the expansion to first order in ui 2 . 



eBHc^ Z" 00 dps K 



where L — 2m 2 c 4 + eBhc(n + n'))F^, + eBinn'Y^G^,. For /3 _1 — > 0, (i.e., for all practical purposes in most 
laboratories) the positron contribution vanishes in (|10I12|) . and we omit it in what follows). For propagation parallel 
to B, kj_ = and the argument i^TI/ vanish. In the approximation assumed above (u> 2 H 2 <C 2e-Bftc)these expressions 
are actually valid also for propagation perpendicular to B. This makes J (n e dp3) / E q the dominant term in the integrals 
in F m , S m . 

c)We return to the Faraday rotation problem. We have 

„ a mc X 3 eBhcN e . . 

F= -2i fa (13) 

For propagation parallel to B we have, 

T m ~ 10" 45 iV e , R = F ~ lO" 48 ^ 
By dividing R by the factor |k| ~ 10~ 12 erg we get 

Aw = R/k ~ 10~ 36 7V e 

Thus, e'/uj = Alu/uj ~ lCT 24 7V e . A density of 10 6 electrons/cm 3 , would give a Faraday rotation frequency whose 
order of magnitude is comparable with PVLAS rotation. 
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